(N 



(N 



O^' 



> 

m 

q 

en 

o 

(N 



% 



Adiabatic approximation in the ultrastrong-coupling regime of a 
system consisting of an oscillator and two qubits 

Ping Yang, Zhi-Ming Zhang^ 

Laboratory of Nanophotonic Functional Materials and Devices, 

SIPSE and LQIT, South China Normal University, Guangzhou 510006, China 

(Dated: February 24, 2013) 
O. 

<^ ■ Abstract 

We present a system composed of two flux qubits and a transmission-line resonator. Instead 
of using the rotating wave approximation (RWA), we analyse the system by the adiabatical ap- 
proximation methods under two opposite extreme conditions. Basic properties of the system are 

Q-c calculated and compared under these two different conditions. Energy-level spectrum of the system 

-(— > \ 

C . in the adiabatical displaced oscillator basis is shown, and the theoretical result is compared with 

^J the numerical solution. 
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I. INTRODUCTION 

Recently, more and more attention has been paid to superconducting devices [l|-[6| to 
build quantum systems, as the tunability of system parameters, which is one of the most 
exciting advantages of superconducting circuit QED over natural-particle-based cavity QED, 
makes such devices more likely to be successful in quantum information processing [7|- 
y]. In fact, some remarkable progresses have^been made in recent years, such as Fock 



states preparation in superconducting devices 10|, single-photon router in the microwave 
regime ll|, and high-fidelity readout in circuit QED [l2|. Superconducting flux qubit that 
consists of superconducting loops and Josephson junctions can be viewed as a two-level 
system when the parameters satisfy the condition called degeneracy point [13|. 

The paper is organized as follows: In Sec. II we introduce the system and the Hamilto- 
nian. In Sec. Ill we discuss the properties of the system by adiabatic approximation under 
two opposite extreme conditions. Energy-level spectrum of the system in the adiabatical 
displaced oscillator basis is shown, and the theoretical result is compared with the numerical 
solution. Sec. IV is the conclusion. 



II. SYSTEM AND HAMILTONIAN 

The system that we consider consists of a harmonic oscillator and two three- Josephson- 
junction qubits (The smallest junction has been replaced by an additional loop [14]) which 
are coupled to the oscillator. Hamiltonian of such qubit is given by Hq^ in Eq. ([2]). 

The oscillator here is a microwave transmission-line resonator, and the qubits which 
are superconducting flux qubits are fabricated so that the loops are closed by the center 
conductor 15|. The schematic diagram of the structure is shown in Fig. [D^a), and it is 
worth noting that the distance between the two qubits is sufficient large such that the 
interaction between them can be ignored. The schematic graph of our system is illustrated 
in Fig. [U^b), in which two qubits are coupled to a harmonic oscillator. It is worth noting 
that the qubits are considered to be identical, which means that the parameters A, e, Eq, 
and the coupling strength g for both of the qubits are of the same value. 

The Hamiltonian corresponding to our system is 
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FIG. 1. (Color online) (a) Schematic diagram of the structure. The two light blue squares are 
improved three-junction flux qubits fabricated to the center conductor, (b) Schematic graph of 
the system. Two identical qubits (i.e. parameters A, e, energy-level splitting Eq and coupling 
strength g for both qubits are of the same value) viewed as two- level system with ground state l^) 
and excited state |e), are coupled to a harmonic oscillator whose characteristic frequency is ujq. 
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Hnt = gX[(T^^ +0-^2J> 

where A is the energy gap tuned by the flux in the additional loop, and e is the bias tuned 
by magnetic flux in the main loop. Energy-level splitting Eq = a/AM-~?, and it will be 
useful for the following analysis to deflne an angle 9 by tan^ = e/A. p and x are the 



moinentum operator and position operator of the harmonic oscillator. (Xj. and dz^ are the 
Pauli operators of the ith qubit. 

We express the Hamiltonian of the oscillator and interaction using a^ and a (creation and 
annihilation operator of the oscillator), so Hog and Hint are rewritten as 

Hos = fyjjQO^a H — hw^^ 

Hint = X{a^ + a)(o-^i +0-^2), (3) 

There is no analytic solution to Eq. ([1]) so far, but some approximation methods have 



ITJI). In the next section, 



been discussed for one-qubit case without RWA (such as |16| . 

we present two adiabatic approximations that can be used to describe our system under 

different regimes of parameters. 

III. ADIABATIC APPROXIMATIONS UNDER TWO OPPOSITE EXTREME 
CONDITIONS 

A. Adiabatic approximation in the displaced oscillator basis 

In this section we consider the case in which Itwq is far greater than Eq. In this case, 
using adiabatic approximation, one can consider that each of the qubits has a well-defined 



value of cr^, i.e., a^^ = ±1, and a^^ = ±1 [l8|. When a^j = a^j = dzl, eigenstates of 



the system can be written as \iIj±, ±) = \iIj±) (g) |±), where|'?/'±) stand for eigenstates of the 
oscillator, |+) and |— ) stand for the qubits' eigenstates |ei,e2) and \gi,g2), respectively. 
When a^i = —0^2 = ±1, eigenstates of the system can be written as IV'ojO) = l^/^o) ® |0), 
where |0) stands for \ei,g2) and \gi,e2)- 

In the case a^^ = 0^2 = ±1 which means the states of qubits are |ei,e2) or \gi,g2), the 
effective Hamiltonian of the oscillator [19] 

Hos,eff\a.,=a.2=^l = huJoo) & ± 2A(a + o)). (4) 

Assuming A and ujq are all real, eigenstates of this Hamiltonian are as follows which can 
be viewed as displaced Fock states: 

1^^) = e^(2A/r«^«)("^-")|n) = l^±), n = 0,1,2,..., (5) 



The eigenenergies are given by 



E' = nhujQ — A\^ /hujQ 



72 = 0,1,2, 



(6) 



In the case a^-^ = —Oz-i = ±1 which means states of qubits are \e\.,g2) or \g\,e2), the 



effective Hamiltonian of the oscillator is the same as the usual harmonic oscillator 



^os,e//U,,=-<7.2=±i = ^o(o^a + 1/2). 



(7) 



So the eigenstates are given by Fock states 



iV'o) = l'^) = I'^o), 



n = 0,1,2, 



(8) 



and eigenenergies are given by 



E:„ = hijjQ{n\ 1/2) 



no 



n = 0,1,2,.... 



(9) 



The eigenstates and eigenenergies obtained by Eqs. 



and (El) build the main 



results of displaced oscillator basis, which will be used throughout the following analysis. 
The potentials corresponding to above results are harmonic oscillator potentials which are 
illustrated by Fig. |2l 

Because the eigenstates of the oscillator depend on the states of qubits now, there will be 
some new properties of the oscillator's states. States of the same potential well (presented 
in Fig. 12]) still maintain the usual orthonormality, that is (m+|n+) = bmn-, ("^ol'^o) = ^mn^ 
{m_\n_) = 5mni while states of different potential wells do not, the displacement operator 
cause a displacement in x, so they are no longer orthogonal to each other. The overlaps 
between displaced oscillator basis of different wells are given by 
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FIG. 2. (Color online) Schematic diagram of the displaced oscillator basis. The horizontal axis 



2muio 



. All three wells maintain the same harmonic character, and usual eigenstates as 



well. The equilibrium position of the left (or the right) well is shifted by a specific constant. The 
shift direction is to the left (or right) when the qubits are in |+) = |ei,e2) (or |— ) = \gi,g2))- 
The middle potential well which is double degenerate corresponds to non-displaced case in which 
the states of the two qubits are opposite, i.e., |0), and the equilibrium position is higher than the 
others. Eigenstates which have the same value of n in the left are degenerate with the right well 
in energy. 



where L^ are the associated Laguerre polynomials. It is worth noting that {rrioln^) = 
(— l)™~"(r7i_|no), and (?7i+|?2o) = (— l)™'~"(T7io|^+), which are useful identities in the later 
calculation. 



Having obtained the eigenstates of the displaced oscillator together with their properties, 
we now focus on the qubits. Taking any specific value of n, one can build an effective 
Hamiltonian of the qubits for that value of n. Due to the fact that for each value of n 
there are four qubits' states, namely |ei,e2) , \ei,g2) , |5'i,e2) , and \gi,g2), the matrix of 
the effective Hamiltonian of qubits is a 4 x 4 matrix in the space defined by |n+, 61,62) , 
\no,ei,g2) , |no,5'i,62) , and \n^,gi,g2). Based on the overlaps between displaced oscillator 
basis obtained previously, we can calculate the elements of this matrix immediately, and the 



matrix is given by 
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q,eff 



/ -e -f(n+|no) 

f(no|n+) 
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(13) 



V -f (ri_|no) -f (n_|no) e / 

where (n+|no), (no|n+), (n_|no) and (no|n_) are calculated to be of the same value as 

Eigenenergies of this effective Hamiltonian are given by 



E'' 



= 0, 



Ko 



(14) 



n 



0,1,2, 



and eigenstates are given by 

|^n_)= (l + 2Ql + 2e„ Vl + 92) \n+, ei, 62) + \n., 91,92) 

+ (Qn + VTTof^ 1^0,61,^2) + (e„ + v/TTel) 1^0,^1,62), 

l^noi) = -|^+, 61,62) + \n_,gi,g2) + 2e„|no, 61,5-2), 

iV'noj) = -1^0,61,52) + 1^0,51,62), 

1^^^+)= (l + 20,^-2e„v^rTe|)|n+,6i,62) + |n_,5i,52) 

+ (^Qn - a/1 + 0^j |no, 61, 52) + (^©n - a/1 + 02 j |no, 51, 62), 
where Or, = 



(15) 



tan 6 



E;^ + E;'^ and E„„ 



e;, + e;V We 



g-2(x/a^o)^L„(4A2/?!i^o)' 

The energies of the system are drawn as -E„± 
present the energy-level spectrum as a function of X/huo in Fig. [3l and these four pictures 
are different in 6. It is obvious that no matter how we set the regime of the parameters, 
Eng remain constant for any specific value of n, so we do not present these horizontal lines 
in the figures. All of the four diagrams are under the condition that huo/Eq = 4 to ensure 
the adiabatical approximation. When A = 0, energy levels of the system are simplified to 
nhuj + Eg which are different from the energy-levels of the usual harmonic oscillator by a 
constant, so they are equally spaced. As A increases, behaviors are evidently different for 
different value of 6. In the case 6 = (i.e. qubits are in the degenerate point), when A 
increases, some slight avoided crossing emerges, and in the limit of large A the levels which 



have the same value of n form pairs. In the case 9 = 7r/6, avoided crossing vanishes and 
sphtting of energy-level pairs with the same value of n turns up in the large A limit. The 
splitting of energy pairs enhances when 6 = tt/4. In the case 6 = vr/S, there is no obvious 
energy pairs any more, and space between each energy levels maintain the same in the 
duration of increasing A. 
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FIG. 3. (Color online) The energy-level spectrum diagram of several lower levels as a function 
of {X/fkvo). hujQ/Eq = 4. The different color denotes different states of the oscillator (black for 
n = 0, red for n = 1 , blue for n = 2 , and green for n = 3). The solid lines stand for qubits in 
state 151,(72)1 and the dashed lines stand for |ei,e2). (a) 6* = 0, in the limit of large A the levels 
which have the same value of n form pairs, (b) 9 = vr/6, avoided crossing vanishes and splitting of 
energy-level pairs with the same value of n turns up in the large A limit, (c) = tt/4:, splitting of 
energy pairs enhances, (d) 9 = it/3, there is no obvious energy pairs any more, and space between 
each energy levels maintains the same in the duration of increasing A. 

We show the comparison between the displaced oscillator adiabatical approximation 
method which is given by Eq. f lT^ and the numerical solution for several lower levels 
in Fig. lU They fit well in almost all regime of A. 

When the qubits are in the degeneracy point, i.e., e = 0, the states of the system core- 
sponding to the eneigys En_, En^ (double degenerate) and En^ are superposition states of 
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FIG. 4. (Color online) The comparison between the displaced oscillator adiabatical approximation 
and the numerical solution for the lowest energy level. They fit well in almost all regime of A. 

the displaced basis, which are given by 

|^n_) = l/4|n+, 61,62) + l/4|rao, 61,5-2) + 1/41^0,51,62) + l/4|n_,5i,52), 
IV'noi) = -l/V2|n+,6i,62) + 1/V2|r2_,5i,52), 
IV'noa) = -I/V^I'^O, 61,52) + 1/721^0,51,62), 

|V^n+) = l/4|n+, 61, 62) - l/4|no, 61, 52) - l/4|no, 51, 62) + l/4|n_, 51, 52). 



(16) 



It is interesting that no matter how we set the regime of the parameters, \ipno ) always 
has the form mentioned in Eqs. (IT5l)-(IT6l). It can be rewritten as \il)no ) = l/v^(— |6i,52) + 
|5'i,62)) ® 1^0), which means that oscillator is decoupled to the states of the qubits. So 
the two-qubit maximally entangled states can be obtained by detecting the state of the 
oscillator. 



B. Adiabatic approximation for the case of high-frequency qubits 

In this section we consider the case that Eq is far larger than huj^. In this case, similar 
to what we have discussed above, usin g ad iabatic approximation one can consider that the 
oscillator has a well-defined value of x p^] . Thus the effective Hamiltonian of the qubits is 



19| 



Hq,eff\x = -7r(o-xi + 0-^2) - T^i^zi + cr,,) + gx{a,, + ct^J. (17) 



The eigenenergies of Eq. (TT7|) are given by 

E,^ = ±y/A^ + i2gx-ey, 

(18) 
Eqo — 0; (double degenerate) 

Unlike above results that energies of high-frequency oscillator are independent of the 
qubits' states, the energies of high-frequency qubits now have a dependence on the position 
X of the oscillator. Thus the effective potential of the oscillator has the form 

K.± = Irnuy ± VA2 + {2gx - e)2, 

I (19) 

where Vos_, Ks^ and Voso corespond to the qubits in states \gi,g2), 1^1,62), and \ei,g2) (or 
15-1,62)), respectively. 

Due to the correction terms in Vos^, when the qubits are in the same state, the effective 
potential of the oscillator is no longer harmonic, while it is the usual harmonic when the 
qubits are in opposite states. However , we can sitll obtain an approximate solution similarly 



to Ref. [l9| 



A renormalized frequency uo is obtained, which is given by 



~2 _ 2 



(20) 



thus the relevant approximate effective potential is given by 



^ 2 ( '^^9 Y r. 

'^-"2'""°4'^^^^^C^J '■ (21) 

It is interesting that the oscillator's frequency now has a dependence on the qubits' states, 
it increases when both of the qubits are in ground states, decreases when both of the qubits 
are in excited states, and stays unchanged when the qubits are in opposite states. 

For Vos+, i.e., the qubits are both in excited states, the approximate effective potential is 
a displaced harmonic potential with an increased frequency, so the states of the oscillator 

10 



are displaced Fock states discussed in Sec. III. A. For Vosq, i-e., the qubits are in opposite 
states, the effective potential is a usual harmonic potential, thus the oscillator's states are 
Fock states. 

It is worthy to note that for Vos_, the renornialized frequency a)o_ turns into imaginary 
when Uq < Ag'^/mEg (the stationary point), which means that under the condition 

^ < 1 (22) 

the system becomes unstable, and the approximate effective potential is not applicative any 
more as the effective potential becomes double- well. In the case of degeneracy point (e = 0), 
by differentiating Vos_ given in Eq. flT9|) . one can obtain that the two minimal points of the 



double- well potential are located at ±xo, where Xq = \/ ^g^ jrrP-ijj^ — A"^ /Ag"^. The energy 
barrier height between the two well is determined by g"^ /itllOq. The effect of introducing a 
finite value of e is that two wells are no longer symmetrical and the location of energy barrier 
is shifted to the left or right. 



IV. CONCLUSION 

In this paper we present a system that consists of two flux qubits coupled strongly to a 
transmission line oscillator, and obtain eigenenergies of the system and the properties of its 
eigenstates. Adiabatic approximation methods under two opposite extreme conditions, i.e., 
the adiabatic approximation in the displaced oscillator basis and the adiabatic approxima- 
tion of the high-frequency qubits, which can be used to analyse our system are compared. 
Although they start by different assumptions, it has been proved that they are both valid 
in the ultrastrong coupling regime. There are some differences between these two approxi- 
mations. It is notable that unlike the later approximation, there is no stationary point that 
turns harmonic potential into double-well potential in the former. 
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